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Abstract 

For an ordered array of critical volatile wetting 
droplets the formation of a super lattice by an 
Ostwald-ripening like competition process is consid- 
ered. The underlying diffusion problem is treated 
within a quasistatic approximation and to first or- 
der in the inverse droplets distance. The approach 
is rather general but a square lattice and a triangu- 
lar lattice are studied explicitly. Dispersion relations 
for the super-lattice growth of these arrays are calcu- 
lated. 

In a recent experiment C. Schafle et al |^ nucleated 
diethylene glycol droplets out of a supersaturated va- 
por phase on a hexagonal array of hydrophilic (alka- 
nethiol) patches on a hydrophobic substrate. The 
pressure of the system was reduced afterwards so that 
the droplets started to shrink by evaporation. Dur- 
ing this process the hexagonal lattice split into two 
triangular super lattices in a sense that the droplets 
on one super lattice shrank faster than those on the 
other super lattice. The super lattice of faster shrink- 
ing droplets disappeared, finally. The remaining tri- 
angular droplet lattice appeared to be stable against 
further development of super structures. The appear- 
ance of the triangular super lattice from the hexago- 
nal lattice is interpreted as an Ostwald-ripening type 
process. The system reduces its free energy by hav- 
ing few large droplets instead of many small [|[ ||, ^ . 
Similar observations where made by the same authors 
101 as well as by Lacasta et al ^ simulating Cahn- 
Hilliard dynamics of droplet arrays in a concentration 
field. The study of super lattice growth on an array of 



wetting droplets is not only of fundamental interest 
as an example of pattern formation and the dynamics 
of wetting on structured surfaces, but it is also im- 
portant to learn about the features of volatile droplet 
arrays since they can be used to build e.g. chemical 
sensors (see ||, 01)- 

The scope of this paper is the analytic discus- 
sion of the diffusional growth properties of a set of 
wetting droplets sitting on hydrophilic patches on a 
hydrophobic substrate. We are considering circular 
patches of equal radius a. Square and triangular lat- 
tice configurations of the patches will be studied ex- 
plicitly. It is assumed that there is a wetting droplet 
sitting in the center of each patch. The different pos- 
sible equilibrium configurations and phase transitions 
in such a system were discussed in great detail by 
Lenz and Lipowsky only recently [|[ § . 

To discuss the diffusional growth of a wetting 
droplet out of the surrounding gas phase, one has 
to specify the boundary conditions of the diffusion 
equation on the substrate and on the droplets surface. 
Since there is no diffusion flux into the substrate, the 
normal derivative of the concentration field vanishes 
there. The concentration Cs in the gas phase on top 
of the droplet is given by the Gibbs-Thomson rela- 
tion. It is a constant along the droplets surface and 
is just a function of the radius of curvature R [0, ^ 



c,{R) = co ( 1 + ^ 



(1) 



where cq is the concentration above a flat condensate 
and A is a capillary length. 
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We assume that the droplet is growing slowly due 
to its high density compared to the surrounding gas 
phase. The diffusion equation, which describes the 
concentration around the droplet, can be approxi- 
mated by the Laplace-equation Ac = 0, therefore. 
The Neumann condition of vanishing normal deriva- 
tive on the substrate suggests to mirror the system on 
the substrate to fulfill the condition implicitly. This 
way the semi-infinite system becomes an infinite sys- 
tem formally and the spherical-cap shaped wetting 
droplet becomes a symmetric lens that provides a 
uniform (Dirichlet) boundary condition to the out- 
side concentration field. Since this field is described 
by a Laplace equation, one can use an electrostatic 
analogy where the electrical potential corresponds to 
the concentration field. The mirrored droplet, i.e. 
the symmetric lens, is a conductor in this picture 
iQ, The volume growth Cl of an individual 

droplet corresponds to the charge of the analogous 
conductor because the diffusion flux density on the 
droplets surface is given by the normal derivative of 
the concentration field and the charge density of the 
conductor is given by the electric field on top of the 
conductor which is the derivative of the potential. 
This argument leads to the growth equation 

n = c{R,e)-{coo-csiR)) (2) 

eventually |l^, where C{R, 9) is the electric capac- 
ity of the symmetric lens and Coo is the systems con- 
centration far away from the droplet. The capacity 
C{R,9) depends linearly on the radius of curvature 
R and is a function of the wetting droplets contact 
angle 6. It can be found in the literature |ll|, [l^ . 

Since the surface tension between the droplet and 
the hydrophilic patch is smaller than the one between 
the droplet and the hydrophobic substrate, the con- 
tact (Youngs) angles are different in these regions. A 
small droplet, which completely sits on a hydrophilic 
patch, has a small contact angle di . If the droplet is 
too big to fit onto the patch it reaches out on the hy- 
drophobic region where its contact angle 02 is bigger. 
Consequently one can distinguish three different vol- 
ume regimes for the droplets. Regime (1), where the 
contact line is on the hydrophilic patch completely, 
regime (2), where the contact line is fixed on the 



border between the hydrophilic and the hydropho- 
bic material and regime (3), where the contact line is 
completely in the hydrophobic region (see fig. |l|). 

In regime (2) droplet growth takes place with a 
fixed contact line but a changing contact angle 9. If 
9 < 7r/2, this leads to the unusual case that the radius 
of curvature R decreases with an increasing droplet 
volume. This implies, that in this case a droplet in 
an environment Coo = Cs(i?) is not critical but sta- 
ble because its Gibbs-Thomson boundary concentra- 
tion increases with increasing volume unlike the usual 
case where the Gibbs-Thomson boundary condition 
decreases with increasing droplet size. This stability 
allows to grow droplets of identical size on an array 
of equally sized patches in the regime (2), which is 
important for experimental purposes. 

Now there is not only one droplet in the system but 
the collective behavior of a system of droplets shall 
be studied. The boundary condition on top of every 
single droplet is given by ([T]), so that we can use the 
full electrostatic analogy again. We have a system 
of conductors where the relation between the charges 
and the potentials is given by a capacitance matrix 
which now yields 

J7j = ^ dj ■ (coo - Cs{Rj)) (3) 
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for the growth rates of the individual droplets gen- 
eralizing eq. (^). Here the indexes i,j count the in- 
dividual droplets, fli and Ri are the volume and the 
radius of curvature of droplet i and the matrix Cij 
is the abovementioned capacitance matrix of the sys- 
tem. 

In the following we will assume that the lattice 
constant d of the droplet array is large compared to 
the droplet sizes {Ri <C d). Then, the inverse of the 
capacitance matrix can be written 

1 _ r C-i(i?, e) + (1/^2) for i = j 
'1 \ for i^j 

where Ci {R, 9) is the capacity of the droplet i alone. 
Rij is the distance between the centers of the droplets 
i and j. Eqs. (0), (|) and (|) lead to 
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to linear order in inverse droplet distances. 

If we now consider a system where each droplet sat- 
isfies Ri = Rc ■— (Aco)/(coo — Co) (i.e. it is critical if 
its radius increases with increasing volume or stable if 
its radius decreases with increasing volume), Ri = 
is a (in most cases unstable) solution of the growth 
equation. To study the onset of the spontaneous de- 
velopment of a super lattice on the regular lattice of 
wetting droplets, we introduce 6Ri • e"* :— Ri ~ Rc 
where the 5Ri describe an eigenmode of the linearized 
version of eq. (^ . The growth exponent reads 



2VL' 



1 SRj 
J^i B^j 5R, 



(6) 



Here Cc is the capacity of an individual critical 
droplet alone and Vl' := dfl{R,e{R))/dR\R^R^ de- 
scribes the variation of the droplet volume when its 
radius R is changed obeying the physical boundary 
conditions on the contact angle, i.e. holding the con- 
tact angle constant if we are in the regime 1 or 3 of 
fig. |l| and pinning the contact line (varying the con- 
tact angle) if we are in the regime 2. As mentioned 
earlier, O' is negative in regime 2 if 9 < tt/2. 

The sum E := d • J2j=it {^1 Rij){5Rj / ^Ri) is inde- 
pendent of the choice of i since we look at eigenmodes 
of eq. (^). It depends on the type of droplet lattice 
and the mode under consideration. One finds a diver- 
gence at E/d = —Cc. This divergence is an artifact of 
the approximation for the capacitance matrix to first 
order in l/Rij, i.e. neglecting 0{l/cP). The fastest 
growing mode is found at the smallest E. In the fol- 
lowing the growth exponent oj will be evaluated for 
the different modes of an infinite square lattice and 
an infinite triangular lattice. 

Due to the symmetry of the square lattice the 
eigenmodes of the growth equation are just Fourier 
modes. In fig. |2| the sum E is plotted as a function 
of m and n, where m, n are the wavenumbers along 
the two axis of the square lattice. The minimum 
is found at E(7r,7r) which means that the preferred 
growth mode leads to a checkerboard structure. The 
two corresponding square super lattices on the orig- 
inal square lattice carry uniformly shrinking or uni- 
formly growing droplets, respectively. This behavior 
is optimal in the sense that every growing droplet is 



surrounded only by shrinking droplets and the other 
way around leading to a fast exchange of matter. 

In case of a triangular lattice things are not as obvi- 
ous since the geometry does not allow every droplet to 
be surrounded by droplets of a different growth sign 
without frustration. If we introduce two non rectan- 
gular axis along two of the primary axis of the trian- 
gular lattice we have again translational symmetry 
along these axis and can Fourier analyze the system. 
We get a sum E for the corresponding wavenumbers 
as shown in fig. ^. The dispersion relation has a three- 
fold degenerate minimum at (7r,7r), (0,7r) and (tt, 0) 
reflecting the special symmetry of the triangular lat- 
tice as shown in fig. |4[ These fastest eigenmodes can 
be characterized by the property that each droplet 
has four nearest neighbors with opposite growth signs 
and two nearest neighbors with the same growth sign. 
In this respect the growth conditions in the triangu- 
lar lattice are inferior to those in the square lattice. 
This is reflected by the fact that the minimum E of 
the square lattice is lower than the one of the trian- 
gular lattice, i.e. the super lattice formation is faster 
in the square lattice case.. 

C. Schafle et al. found in their experiment Q that 
the triangular droplet lattice that occurred as a super 
lattice of the initial hexagonal lattice did not show 
any further instability, whereas the calculations pro- 
vided in this paper display an instability as sketched 
in fig. H Unlike the subject of this paper the ex- 
periment does not analyze critical droplets so that in 
addition to oj there is a timescale related to the over- 
all evaporation of droplets. It seems reasonable to 
expect that the instability of the triangular droplet 
lattice could be seen experimentally by stretching the 
timescale provided by the overall evaporation, i.e. by 
going closer to the critical concentration. 
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Figure 2: The sum S as a function of the wavenum- 
bers m and n for a square lattice. 




[12] C. Snow, J. of Research NBS, 43, 377 (1949) 



Figure 3: The sum E as a function of the wavenum- 
bers m and n for a triangular lattice. 




Figure 1: Three qualitatively different droplet sizes 
can be distinguished: (1) A droplet which sits on the 
hydrophilic patch completely; (2) the contact line of 
the droplet sits on the edge of the hydrophilic region; 
(3) the contact line is on the hydrophobic substrate 
completely. 
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Figure 4: The fastest growing modes of a triangular 
lattice. 
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